Adhesion Induced DNA Naturation 
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DNA adsorption and naturation is modeled via two interacting flexible homopolymers coupled 
to a solid surface. DNA denatures if the entropy gain for unbinding the two strands overcomes 
the loss of binding energy. When adsorbed to a surface, the entropy gain is smaller than in the 
bulk, leading to a stronger binding and, upon neglecting self-avoidance, absence of a denatured 
phase. Now consider conditions where the binding potentials are too weak for naturation, and the 
surface potential too weak to adsorb single strands. In a variational approach it is shown that their 
combined action may lead to a naturated adsorbed phase. Conditions for the absence of naturation 
and adsorption are derived too. The phase diagram is constructed qualitatively. 

PACS numbers: 82.35.Gh, 82.39.Pj, 05.90.+m 



The structure of biopolymers is the key for understand- 
ing their biological functioning 0, ■ A known compo- 
nent of this structure is the double-strand — two polymers 
lined up by hydrogen bonds — realized in DNA and re- 
sponsible for stability of the genetic information. 

We shall study the adsorption and surface (adhesion) 
induced naturation of a double-strand polymer with re- 
solved motion of the separate strands. The adsorption of 
DNA onto surfaces is of interest for a number of applica- 
tions such as design of substrates as carriers of genetically 
engineered vaccines, surface patterning, histone-DNA in- 
teractions, and electrophoresis 0. DNA adsorption can 
be realized by different mechanisms 0: covalent bond- 
ing, electrostatic attraction to a positively charged sur- 
face jjj, hydrogen bonding, and hydrophobic forces. 

When the two strands are tightly connected to each 
other, the double-strand is described via a single semi- 
flexible chain d, H, • Indeed, for the double-stranded 
DNA the corresponding persistence length l p is about 100 
nm 8], or 300 base-pairs. However, this large value of l p 
is more a result of "construction rigidity" , because for a 
single DNA strand l p is much smaller, about 1 nm (2-3 
base-pairs) 0, while the length of the strand is 10-100 
fim. The single DNA strand can thus be modeled as a 
flexible chain [jj. Adsorption of a (semi)flexible chain is 
by now a classical problem in polymer physics 0- 0- 01 ■ 

If the motion of the single strands is resolved, i.e., when 
the inter-strand hydrogen bonds are relatively weak, the 
polymer becomes a complex system with different, mutu- 
ally balancing features. A realistic model of DNA should 
include stacking energy between two base pairs, the he- 
lical structure, self-avoidance, heterogeneity, etc. Such 
models do not seem to exist; there are, however, various 
models with different degrees of sophistication c apturing 
at least some aspects of the problem palalrj|^ ll(All^ll2| . 

Our model deals only with most basic features of the 
problem. We consider two flexible Gaussian chains in- 
teracting with each other and with a plane solid surface. 



Though neglecting some features of DNA (see below), 
our model predicts two mechanisms of adhesion induced 
naturation. This model without surface was mentioned 
in [j| and studied in in the context of DNA denatu- 
ration (melting); see also f° r recent results. 

Consider 2N coupled monomers with radius-vectors 
fi|fe and ?2\k (k — 1, N) and potential energy 



fc=i 
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2 \r<x\k ~ r a \k-i)" + V(f a |fc) 



where fk = fW — f^ife, hi is the inter-strand potential 
and V is the surface-monomer potential. The harmonic 
interaction with stiffness I is responsible for the linear 
structure of the polymers. The system is embedded in 
an equilibrium thermal bath at temperature T = 1//3 
(fee = 1). The quadratic kinetic energy is irrelevant, 
since it factorizes from the partition function and does 
not influence the monomer density. It is known that in 
the thermodynamical limit N 1 the free energy of a 
single flexible chain (without self-interactions) is deter- 
mined from an effective Schrodinger equation 00. The 
considered two-strand situation is a direct generalization 
of this, with the two-particle equation being = Ety , 

±d 2 + V(z a )] + U(r), d a = d fa , (1) 

Here V = l(3 2 V, U = lf3 2 U, while r xa = (^,2^1,2,^1,2) 
are the position vectors of two quantum particles repre- 
senting the strands and r = r\ — ri is the mutual position. 
We shall extensively use this quantum language. If there 
is a gap between the lowest two eigenvalues of H, the 
ground state wave-function ^ determines the monomer 
density as n(f{, rj) = *if 2 (ri, r 2 ). The eigenvalue E is the 
energy of the quantum pair, related to the free energy F 
of the polymer as E = lj3 2 F/(2N). 

The surface is described by an infinite potential wall, 



H 



*(ri,f 2 ) = 0, if z\ < 0, or z 2 < 0. 



(2) 
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Both V(z) and U(r) are attractive, V < 0, U < 0, and 
short-ranged: / °° dzV(z) and J °° drr 2 U(r) are finite. 
When J7 = 0, the Hamiltonian H reduces to two uncou- 
pled strands (or one-dimensional particles), each in the 
potential V{z). If V{z) is shallow enough, no bound (neg- 
ative energy) state exists, while the second-order binding 
transitions corresponds to adsorption of a single flexible 
polymer Q. The physical order-parameter for this tran- 
sition is l/(z 2 ), which is finite (zero) in the adsorbed 
(desorbed) state. Denote by p the dimensionless cou- 
pling constant of V = pV such that (for U = 0) the 
adsorption threshold is p c = 1. Note that the adsorption 
of a single strand DNA is a part of the renaturation via 
hybridization , a known method of genetic systematics. 

Analogously, switching off both V(z) and the wall, we 
shall get a three-dimensional central-symmetric motion 
in the potential U (r) which again is not bound if U is 
shallow. This second-order unbinding transition with the 
order parameter l/(r 2 } corresponds to thermal denatu- 
ration (strand separation) of the double-strand polymer 
Writing likewise U = XU, A is the dimensionless 
naturation strength. We take the naturation threshold 
in the bulk to be A c = 1. When the wall is included, i.e. 
© is imposed, the strands loose in the adsorbed phase 
part of the entropy needed to denaturate. 

As seen from Eq.£Q), in the (x,y) plane H is both ro- 
tationally and translationally symmetric. Thus, for the 
low energy physics the only relevant variable coming from 
that plane is p — y/ (xi — x 2 ) 2 + (yi — y 2 ) 2 - The relevant 
wavefunction is thus ip(zi, z 2 , p) and the reduced Hamil- 
tonian H r reads 



H r = H 1 +H 2 - p-'dppdp + U (v / P 2 +^ 2 ) , (3) 

(4) 



H ia ^~-d 2 Zi2 +V{z ia ) 



z = Zl - z 2 . 



There is no hope for an exact solution of this model, since 
it belongs to the class of three-body problems (the role 
of the third particle is being played by the surface; see 
also below). We shall estimate the lowest energy Eq of 
H r via a variational principle: Eq < (x\H\x), where x 
is a normalized trial function. Assume that V{z) (with- 
out U) has a lowest (negative) energy state with normal- 
ized wave function </>(z): H\ t 2^(^1,2) — -EV<?K z i,2), with 
0(^1,2 < 0) = 0. We then make the Ansatz 

r- OO 

X(p,zi,z 2 ) =u}(zi,z 2 )£(p), / dpp£ 2 {p) = 1. (5) 

Jo 

where u>(zi,z 2 ) — <p(zi)(j){z 2 ) and where £(p) is normal- 
ized. Varying (x\H — E\x) w.r.t. £(p), we get an effective 
two-dimensional problem: —p 1 d p (pd p £,) + U c s(p)£, = 
(E — 2Ey)£, where by definition E > E , and where 
U c (p) is an effective potential: 

f°° 2 
U c (p)= [/( vV + (z! - z 2 ) 2 )TT dz a cf> 2 (z a ). (6) 

JO J-A Q= 1 



Like U, U c is attractive and short-ranged. It is well 
known that any (no matter how weak) attractive poten- 
tial in two dimensions creates a bound state, though the 
localization length of this state is exponentially large for 
small energies 131. Thus there is a normalizablc function 
such that E < E < 2Ey < 0, i.e., there is an overall 
bound (double-strand and localized next to the surface) 
state provided V(z) creates a bound state. In the present 
model strong enough surface potential prevents denatu- 
ration (melting) of the double-strand. Upon taking into 
account the neglected features (helicity, self-avoidance or 
stacking), the denaturation transition temperature in the 
adsorbed phase is enhanced but (presumably) finite. 

Next to the threshold p, c = 1 of V(z), 2Ey = ~n 2 is 
small and the wave-function is almost flat [13j: <p{z) ~ 
y/2Ke~ KZ . Thus Eq. © implies that U c = O(k), so 
naively it predicts absence of binding at the threshold 
k = of V(z). But let us study the region where V{z) 
alone does not create any bound state, by employing the 
trial function JSJ, where now ^(p) is a given normalized 
function (its form is irrelevant for us), while u(zi,z 2 ) is 
found from the variational equation: 

[H 1 + H 2 + U eS (\z 1 -z 2 \)-E]u;(z 1 ,z 2 )=0, (7) 



U eS (z)= dp P e( P )U(Vp T +^), (8) 
Jo 

Equations JJJ and Q describe two one-dimensional par- 
ticles interacting via short-range potential U c s = \U c g. 
We shall show that Eq. JJJ) predicts an overall bind- 
ing — that is it predicts E < and a localized normal- 
izable wave-function u>(zi, z 2 ) — at the threshold k = of 
the potential V(z). For a large A this is expected, since 
two strongly coupled particles roughly behave as one with 
double mass in a double potential, which decreases the 
threshold value of V{z) by a factor 4, i.e., /i c (oo) = j. 
Next we focus on small A. Though the energy 2Ey is 
nearly zero, one can apply an ordinary perturbation ex- 
pansion in A, since the suitable matrix elements of U e s 
appear to be small as well, grace to the small prefac- 
tor k in the zero-order wave-functions |00) = (f>(zi)(f>(z2). 
Noting that (00|C/ cff |00) = 0(k) for k -> and thus neg- 
ligible, we get at the second order 0] 



E 



dK 



|(00|£7eff|g)| ; 
SK + K 2 



(9) 



where the integration over the collective variable K 
involves all excited wave-functions of the unperturbed 
two-particle system with wave-vector K and energy 
Ek- There are three orthogonal families of these states: 

|ofc) = 4>{zxXk{kz 2 ), |fco) = Ck(ft«iM*2), 1*1*2) = 

Cfci (/ci2i)(fc 2 (k 2 z 2 ), where Ck(kz) are the correspond- 
ing single-particle excited (continuous spectrum) wave- 
function of the potential V(z) with the wave-number k. 
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For k — > we get after some steps the finite result 

n 2 



E 



= -128 



dk 



o 1 + ^ 2 
dfcidfc2 



dv 



-3v 



Co(kv) 



l 

2 



fcf 



^2 



(10) 
n 2 



dve v (o(k 1 v)( {k 2 v) 



where S = J °° duU c ff(u). With (o(kv) = y § sin(£;td we 

get from I|1U|I the numerical value E = —0.45 S 2 . Its order 
of magnitude can be checked from the exact relation: 



(00|Z7 eff |cj) = (E - 2Ev)(00\c 



2E V = —K 



(11) 



Assume that V(z) is close to its threshold (i.e., [i — > 1 + , 
or k — ► 0) and that a weakly-bound state exists with 
energy E = —p 2 . Recall that 1 00) is mainly a constant 
modulo a numerical factor which cancels from . Since 
all potentials in (0) are short-range and thus negligible 
for large z\ and z 2 , we get to(zi, z 2 ) = f(zi, z 2 )e~ plyZl+Z2 \ 
where / — ► 1 for z\, z 2 — > oo. This is put into and 
for /c — > and small p one gets indeed E ~ — S 2 ~ —A 2 . 

For < 1 and A < 1 we are in the situation where nei- 
ther V^(zi) +V(Z2) (the attractive wall alone) nor U alone 
can create a bound state. We thus conclude from (|10llll|) 
that the present approach does predict binding for fj, = 1 
and for sufficiently small A. Since the ground state is 
supposed to be continuous in A and /i, the very fact of 
having a negative energy at fi — 1 and not very large A 
implies that a naturated, adsorbed state will also exist 
for /i c (A) < fi < 1, where neither of the potentials V and 
U alone allows binding. The precise curve /i c (A), where 
the ground state energy is equal to its value at V = 
(adsorption threshold), requires numerical analysis to be 
reported elsewhere . We thus have found an example 
of so called Borromean binding 0], where the involved 
potentials do not produce bound states separately, but 
their cumulative effect does so. We see from (0 that this 
unusual binding is connected to correlations between the 
^-components of the particles and (separately) to corre- 
lations between their x and y components (in p). 

No-binding conditions. To complete the phase dia- 
gram, we employ a method suggested in [l5J. Let us 
introduce a third Active particle with mass M and the 
radius vector r 3 = (£3,2/3,-23). This particle will substi- 
tute the surface: V(z a ) (a = 1,2) becomes V(\z a — z 3 |), 
and condition (J5J is put at z a — z% only. The modified 
problem with Hamiltonian Hm reduces to the original 
one when for M — ► 00, F3 becomes a fixed vector which 
can be chosen at the origin. We have 

Hm = ~2M + U{r) + ^ + V(lZa Z3|)] - 

It is invariant when shifting the radius vectors f a (a = 
1,2,3) over any fixed vector, while for a finite-particle 
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FIG. 1: Schematic phase diagram for the inter-strand cou- 
pling A versus the strand-surface coupling /i. The bold lines 
confine three thermodynamical phases. ND: Naturation and 
desorption. NA: Naturation and adsorption. DD: Desorp- 
tion and denaturation. The critical naturation strength in 
the bulk is A c = 1, for single strand adsorption it is fi c = 1. 
The following subregions are confined by normal lines, a: Do- 
main described by the no-binding condition I16H . b (bounded 
by the bold DD-NA line and two straight segments): Bor- 
romean naturation and adsorption, c: Adsorption and natu- 
ration due to overcritical coupling to the surface. 



quantum system, a symmetry of Hm implies the same 
symmetry for its ground-state wave- function |l3|. Thus, 

* = #(r! - ra, n - f 3 , r 2 - r 3 ), Y] 3 d a ^ = 0.(12) 

— J a— 1 

We now decompose Hm = #o + #12 + #13 + #23, where 

2d 3 + bJ2 2 , ^1 -]T 3 B a , (13) 

^ — ^l — 1 ^ — 'a— 1 

/ -. -, \ 2 

c I d\ — xd 3 



He, 



H 



13 



1 + X 



+ V(\zi-z 3 \), (14) 



Hli = -2dl^—^] +U(\n-f 2 \), (15) 



while H 2 3 = Hi3(l — » 2). The coefficients a, b, c and d 
can be read off directly: a = toM -1 — 2x 2 /(l + 2a;) 2 , b = 



d = 2a;(l-|-2;)/(l + 2ir) 2 



(1 +x)7(l + 2a;) 2 . Since 



c and d are to be employed as effective masses, for the 
free parameter x we consider x > (the regime x < — 1 
is of no help). Note that ($\Ho\&) =0 due to lfT2)l. Let 
us see when (^|i?^f|^) > 0, i.e, when a bound state is 
absent. Changing the variables — e.g., s*i = (1 + x)fi, 
s 3 = (1 + x)n + f 3 for Hi 3 — we see that (*|F 13 |*) and 
(^liJia |^) are positive for, respectively, /i < c(x) and 
A < 2d(x) . There is no binding under these conditions for 
any M including M — > oo which returns to the original 
problem. By varying x from to oo, we cover the domain 



1 



0</i< -(2- A + 2V1- A) 
4 



(16) 



This is the no-binding region a shown on Fig. ^ At 
A = 1 one needs at least /i c = 0.25 to achieve binding, 
while at fi = 1, A = + suffices. 
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As shown on Fig. ^ there are three phases in the 
model: denaturated and desorbed (DD), naturated and 
adsorbed (NA), and naturated and desorbed (ND). 
The denaturated and adsorbed phase is present only for 
U = 0; recall the discussion below The exact bound- 
ary between DD and NA phases has to be convex: if 
(Ai, /ui) and (A2,^a) are in DD — i.e., H(Xi,fj,\) > 
and H{\<i,H2) > — then the entire line joining these 
points is in DD, because A and \i enter H linearly: 
H(e\i + (1 - e)A 2 , e/ii + (1 - e)jU 2 ) > 0, where 1 > e > 0. 

We already mentioned that for large A the adsorption 
of naturated DNA will happen at fx c = | due to doubling 
of the mass and potential. For A approaching A c = 1 
from above, the absorption is weaker and the entropy will 
increase so we expect that the potential U will facilitate 
adsorption less, leading to the ND-NA curve with fi c > 
j, see Fig. 1. At A = 1 it is expected to meet DD-NA 
at a point, that separates the three phases ND, NA and 
DD. The DD-ND line is implied by the fact that in the 
desorbed phase the polymer is far from the wall. 

In conclusion, we modeled the surface adsorption of a 
double-stranded DNA. There are several situations where 
the description of this process via a single chain will not 
be adequate, and the two-chain modeling is needed. (1) 
The adsorption realized by the same hydrogen-bonding 
mechanism as the naturation [3j . Then in the vicinity of 
the (temperature or pH induced) melting transition also 
the binding to the surface will be small. Thus the motion 
of separate chains will be resolved. (2) The binding to hy- 
drophobic surfaces (e.g., aldehyde-derivate glass) goes via 
partial melting which exposes the hydrophobic core of the 
helix and leads to the DNA-surface attraction Both 
naturation and adsorption are simultaneously weakened 
by increasing the pH (3) For homogeneous DNA 

at normal conditions (pH= 7 and NaCl concentration of 
0.15 M) the thermal melting occurs at temperatures 67 C 
and 110 C for A-T and C-G unbinding, respectively 
This temperature can be radically decreased by increas- 
ing the pH factor, i.e., by decreasing the concentration of 
free protons in the solvent, since the negatively charged 
phosphate groups on each strand are not screened any 
more by protons and strongly repel each other yj . Due to 
the same reason, for the DNA adsorption on a positively 
charged surface [4j, the increase of the pH will increase 
the electrostatic attraction to the surface. Thus weakly 
naturated and weakly adsorbed states may be produced 
by controlling the pH and the surface charge. 

Our model describes DNA strands as flexible ho- 
mopolymers. Many features of real DNA are put aside 
in this way: stacking, self-avoidance, helicity, hetero- 
geneity of the base-sequence, dependence of the surface- 
strand interaction on the denaturation degree (in the nat- 
urated state the bases do not participate in the surface- 
strand interaction), etc. However, the model can still 
be useful in clarifying the mechanisms of DNA adsorp- 
tion/naturation. We saw that a sufficiently strongly at- 



tracting surface prevents thermal denaturation (melting) 
of the double-strand, because effectively the problem re- 
duces to binding in two dimensions, where any attractive 
potential creates a bound state. In practice this implies 
that in the presence of the surface the melting temper- 
ature increases. It was shown that the attractive sur- 
face and inter-strand coupling together can create bound 
states (i.e., double-stranded and adsorbed to the surface), 
even when neither of these interactions alone is capable 
of creating bound states. This unusual type of binding is 
related to finely correlated motion of the strands. Such 
Borromean states, where binding is due to the cumula- 
tive effect of several potentials, first appeared in nuclear 
physics 0, 0] , but their detection in that field is dif- 
ficult. Their observation in polymer physics might be 
easier, since the involved scales are mesoscopic. Our re- 
sults led to the schematic phase diagram which, as we 
argued, may be checked experimentally. 
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